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Abstract 

When the mass difference between the hghtest slepton, the NLSP, and the hghtest neutrahno, 
the LSP, is smaher than the tau mass, the hfetime of the hghtest slepton increases in many orders 
of magnitude with respect to typical lifetimes of other supersymmetric particles. These small mass 
differences are possible in the MSSM and, for instance, they correspond to the coannihilation re- 
gion of the CMSSM for ^ 700 GeV. In a general gravity-mediated MSSM, where the lightest 
supersymmetric particle is the neutralino, the lifetime of the lightest slepton is inversely propor- 
tional to the square of the intergenerational mixing in the slepton mass matrices. Such a long-lived 
slepton would produce a distinctive signature at LHC and a measurement of its lifetime would 
be relatively simple. Therefore, the long-lived slepton scenario offers an excellent opportunity to 
study lepton flavour violation at ATLAS and CMS detectors in the LHC and an improvement of 
the leptonic mass insertion bounds by more than five orders of magnitude would be possible. 
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I. INTRODUCTION 



Cosmological observations have confirmed the existence of the non-baryonic dark matter. 
The observed dark matter relic abundance is QoMh'^ — 0.11 and points to the existence of a 
stable and weakly interacting particle with a mass m = 100 — 1000 GeV. In supersymmetric 
(SUSY) models with conserved R-parity, the Lightest Supersymmetric Particle (LSP), usu- 
ally the lightest neutralino, is stable and it is a perfect candidate for the dark matter. In fact, 
although the dark matter abundance in an arbitrary point of the SUSY parameter space 
is usually too large and only in regions where the LSP abundance is lowered effectively by 
annihilation the observed abundance can be reached. One of the best mechanisms to lower 
the dark matter abundance is the so-called coannihilation [1]. In the coannihilation region, 
the Next to Lightest Supersymmetric Particle (NLSP) has a mass nearly degenerate to the 
LSP. Then the NLSP and the LSP decouple from thermal bath almost simultaneously and 
the LSP annihilates efficiently through collisions with the NLSP. The degeneracy required 
for the coannihilation to occur is generally (5m/mLSp < a few %, where = /wnlsp ~''^lsp- 
Numerically, this mass difference at a fixed tan (3 value ranges from Sm ~ 10 GeV for low 
Ml/2 to practically zero at large M1/2, where M1/2 is the universal gaugino mass. 

Furthermore, it was pointed out in [2, 3, 4] that if staus are the NLSP and its mass 
difference with the LSP is less than the tau mass, they destroy ''Li/'^Be nucleus through the 
internal conversion. With this small mass difference, staus become long-lived, they survive 
until Big-Bang Neucleosynthesis (BBN) starts and form bound states with nucleus. The 
stau- nucleus bound states decay immediately by virtual exchange of the hadronic current. 
In this way, it was shown that the relic abundance of the light elements is lowered effec- 
tively and the discrepancy between the observed value of ^Li/^Be abundance [5, 6] and the 
predicted value from the standard BBN [7] with WMAP data [8, 9] can be solved. Thus, a 
scenario with stau NLSP and neutralino LSP with 6m < rrir could explain the relic abun- 
dance of the light elements as well as the abundance of the dark matter. At this point 
we could ask whether this small mass difference is really possible in the parameter space 
of the Minimal Supersymmetric Standard Model (MSSM) or even the Constrained MSSM 
(CMSSM) consistently with dark matter constraints. As we will show below, 6m < m^- and 
even smaller values are possible in the CMSSM. As shown in [2], if the mass difference, 6m, 
between neutrahno and stau is less than the tau mass, the two-body decay is forbidden and 
the lifetime of stau increases by more than 10 orders of magnitude. 

Strictly speaking the previous discussion is correct only in the framework of the MSSM 
without intergenerational mixing in the slepton sector, as could be, for instance, the CMSSM 
defined at the GUT scale without right-handed neutrinos, or neutrino Yukawa couplings. In 
this case, the NLSP is only a combination of left and right-handed staus. However, there 
is no doubt that, in a general MSSM we naturally expect some degree of intergenerational 
mixing in the sfermion sector. Even starting from a completely universal MSSM at Mqut 
in the presence of neutrino Yukawa couplings, a small mixing between different sleptons is 
generated by renormalization group equations. Moreover, there is no fundamental reason 
to restrict the flavour structures to the Yukawa couplings and keep the soft mass matrices 
universal. In fact, we would expect the same mechanism responsible for the origin of flavour 
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to generate some flavour mixing in the sfermion sector. As an example, we have flavour 
symmetries [10, 11] where we expect similar flavour structures in the Yukawas and the 
scalar mass matrices. Therefore, in general, we expect that the lightest slepton, li, is not 
a pure stau, but there is some mixing with smuon and selectron. In the presence of this 
small intergenerational mixing, even with 6m < m^, other two-body decay channels like 
^1 Xi + ^(/i), are still open [2]. In this case, given that the flavour conserving two- 
body decay channel is closed, the slepton lifetime will have a very good sensitivity to small 
Lepton Flavour Violation (LFV) parameters. The lifetime will be inversely proportional 
to the sfermion mixing until two-body decay width becomes comparable to three or four- 
body decay widths. Thus, it is worthwhile to study the lifetime of the sleptons in small 
dm case to obtain information on LFV parameters. Several papers have studied in the past 
the possibility of measuring LFV at colliders [12, 13, 14, 15, 16, 17, 18, 19]. However, 
in these papers the LFV decays are always subdominant with small branching ratios and 
is never possible to reach the level of sensitivity we reach in our scenario. Only Ref. [16] 
considers LFV in e~^e~ colliders through the slepton production in a long-lived stau scenario 
in a gauge-mediated model and as in the other works they are not sensitive either to the 
presence of small lepton flavour violation that we consider here. 

Long-lived charged-particles are very interesting since they provide a clear experimental 
signature at the LHC [20, 21, 22]. In Ref. [23] they conclude that even if the lifetime of the 
decaying particle is much longer than the size of the detector, some decays always take place 
inside the detector and it is possible to measure lifetimes as long as 10~^ — 10"'^ seconds in 
a particular gauge mediation scenario. 

In this paper, we study LFV in the coannihilation region of the MSSM with neutralino 
LSP and slepton NLSP with a mass difference dm < rrir. The paper is organized as follows. 
In Sec. 2, we study the case of the MSSM without flavour mixing in the slepton sector. 
Here, we show that a sizable part of the coannihilation region at large M1/2 has 6m < m^ 
region satisfying all the experimental constraints and we calculate the lifetime of the NLSP. 
Then, in Sec. 3, we introduce a source of lepton flavour violation in the slepton mass 
matrices and we recalculate the lifetimes in terms of the LFV parameters. In Sec. 4, we 
discuss the expected phenomenology of this scenario at the LHC and how LHC data could 
be used to measure small LFV parameters. Finally in section 5, we present our summary 
and discussion. 



II. LONG-LIVED STAU IN MSSM 

We start our analysis looking for allowed regions in the CMSSM parameter space where 
the mass difference, (5m, between the lightest neutralino mass, m-^o, and the lighter stau mass, 
mfj, is smaller than tau mass, m,-. In this section, we consider the case of a constrained 
MSSM defined at the GUT scale without neutrino Yukawa couplings. Thus, in this case, fi 
is equal to a combination of the right and left-handed stau and no lepton flavour violation 
(LFV) effects are present. The CMSSM is parametrized by 4 parameters and a sign, 

{Ml/2, "^0, ^0, tan/3, sign(yu)}, (1) 
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where the first three parameters are the universal gaugino and scalar masses, and the uni- 
versal trilinear couplings at the GUT scale, tan j3 is the ratio of vacuum expectation values 
of two Higgses and /i is the Higgs mass parameter in the superpotential. In our numerical 
analysis, we use SPheno [24] to obtain the parameters and mass spectrum at electroweak 
scale and use micrOMEGAs [25] to calculate relic dark matter abundance. 

First, we list the experimental constraints that we have used in our analysis. The main 
constraint, apart from direct limits on SUSY masses, comes from the relic density of the 
dark matter abundance as reported by WMAP collaboration in Ref.[9]: 

VtY^y^h^ = 0.1099 ± 0.0062. (2) 

In this paper, we use more conservative range, 

0.08 < ^DM/i^ < 0.14, (3) 

but the results would remain basically unchained with a smaller range. 

A second constraint comes from the anomalous magnetic moment of the muon, a^ = {g — 
2)^/2, that has been measured precisely in the experiments at BNL. The final experimental 
result [26] is 

Qjexp) _ ggg 208.0 ± 6.3) X 10"^°, (4) 
and the standard model prediction is 

q(sm) ^ 1^^^ ggg 5 ± 5 X 10"^°. (5) 

Therefore, the difference between the experiment and the standard model prediction is given 
[27] as 

Aa^ = G^T^P) - a^f*') = (27.5 ± 8.4) x 10"^° (6) 

which corresponds to a difference of 3.3 standard deviations. These results seem to require 
a positive and sizable SUSY contributions to a^, although it is still not conclusive. In the 
following, in all our plots we will show the regions favored by measurement of at a given 
confidence level, but we will not use it as a constraint to exclude the different points. The 
bound on Br(6 — >■ 57) is also taken into account, but as we will see, in the region of small 5m 
and correct dark matter abundance, M1/2 > 700 GeV. Therefore the numerical prediction 
we obtained for Br(6 — > 37) is always very close to the SM predictions. 

In Fig. 1, we show the allowed parameter regions in the M1/2 - fno plane for tan/3 = 
10, 30, 45 and Aq = 600 GeV. The sign of fi is taken to be positive in all figures as 
required by the b s'~f and constraints. In these figures, the solid lines indicate the 
different confidence level regions for the constraint. The yellow (light) narrow band is 
the region with correct dark matter abundance corresponding to the coannihilation region 
and with dm > rrir- The red (dark) band in the coannihilation region represents the points 
with correct dark matter abundance and 6m < m^-. The green regions are not consistent 
with the dark matter abundance, Eq. (3), while the white area is the stau LSP region and 
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FIG. 1: The allowed parameter regions in Mx/2 ~ plane fixing Aq = 600 GeV. tan (3 is varied 
10, 30, 45 shown in each figure. The red (dark) narrow band is consistent region with dark matter 
abundance and 5m < rrir and the yellow (light) narrow band is that with 5m > m-r. The green 
regions are inconsistent with the dark matter abundance, and the white area is stau LSP region 
and excluded. The favored regions of the muon anomalous magnetic moment at la, 2a, 2.5a and 
3a confidence level are indicated by solid lines. 

therefore is excluded. Note that in models where the observed dark matter abundance is 
accounted for other dark matter components besides the neutralino, the small Sm region 
could be extended to the green regions below the yellow band. Notice that choosing a 
positive and large Aq ( < Sijiq), the renormalization group evolved value A{M]v) is reduced 
with respect to the cases of zero or negative Aq. Therefore the left-right entry in the stau 
mass matrix, A — fi tan jS, is also reduced. Then the region of neutralino LSP, for a fixed 
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TABLE I: Approximate mass ranges for points in the coannihilation region with ruf^ — m^o < 1.77 
GeV and within the 2a allowed region for for tan (3 = 30, Aq = 600 GeV. These points correspond 
to mo G (230,250) GeV and M1/2 G (700,800) GeV. u, d and e, z>i denote the first two generation 
up-type, down-type squarks, and sleptons and sneutrinos. 

value of the neutralino mass, fn^o (~ O.4M1/2), corresponds to smaller values of mg for 
positive Aq than for zero or negative Aq. Therefore this implies slightly lighter spectrum 
and larger contributions to a^. Anyway, the main features of the plots are maintained 
for different values of Aq with only this slightly increase in the masses. From these plots 
we can see that the interesting region of correct value of dark matter relic abundance and 
5m < rrir corresponds always to relatively large values of M1/2. This is due to the fact that 
the coannihilation cross section decreases at larger SUSY masses [28] and for sufficiently 
large M1/2 the cross section is too small to produce a correct dark matter abundance, even 
for TUf^ ~ ^xl' 

Comparing the different tan/3 values, we can see that at low tan/3 the long-lived stau 
region corresponds to smaller values of mo ~ (120, 200) GeV, although it can only reach the 

favored region at 3 a. For tan (3 = 30, the long-lived stau region corresponds to similar 
values of M1/2, but the required values of mo are nearly a factor two larger. However, in 
this case the larger value of tan/3 allows this region to reach the favored region at the 
2a level. At tan/3 = 45 the long-lived stau region corresponds to larger values of both M1/2 
and mo and it can only reach the favored region at 2.5 a. This behavior of the allowed 
regions can be understood as follows. Since the SUSY contribution to is proportional to 
tan/3/m|o, it is small when tan/3 = 10. But it is also suppressed at tan/3 > 30 because 
the mass of the neutralino becomes heavier for large tan/3. With fixed Aq, the neutralino 
mass is determined mainly by M1/2 and the increase of the neutralino mass according to the 
increase of tan/3 is seen in Fig. 1. 
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(a) 
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FIG. 2: Feynman diagrams of stau decay: (a) two-body decay, (b) three-body decay, (c) four-body 
decay. I denotes the charged leptons. 



In these figures, we have seen that Sm < nir is possible in the CMSSM at M1/2 > 700 GeV 
and mo > 200 GeV. For these parameters, the resulting SUSY mass spectrum is relatively 
heavy. In Table I, we show the SUSY mass spectrum for the interesting region at tan f3 = 30 
and Aq = 600 GeV. As can be seen in this table, masses of the lightest neutralino and the 
lightest slepton are around 300 GeV. It is important to emphasize that the lightest neutralino 
is mainly bino and the second lightest neutralino is wino with a small admixture of higgsino. 
The second lightest neutralino and the lightest chargino have masses ~ 550 GeV. The 
remaining two neutralino states and the heavier chargino are mainly higgsinos with masses 
~ 800 GeV. The right-handed sleptons are in the range 300-390 GeV and the left-handed 
sleptons have masses of 510-570 GeV. The squarks are in the range 1150-1560 GeV while 
gluinos are heavier than the squarks and about 1540 — 1680 GeV. Finally the lightest Higgs 
has a mass of ~ 116 GeV and the heavier Higgses have masses of 800-900 GeV. With this 
spectrum, the lightest slepton is produced through decays of the heavier neutralinos, the 
charginos and the left-handed squarks in longer decay chains. Neutralinos and charginos are 
produced via decays of squarks and gluinos. Although squarks and gluinos are relatively 
heavy, they are still in kinematical reach and can be produced at LHC experiments. 

The decay rates of the NLSP stau into two, three and four bodies shown in Fig. 2 are 
approximately given by 



J- 2 -body 
Ts-body 
r4-body 



tan^ Ow 
27rm?, 



6m \/ (5m)^ — m 



di^pf^ cos"' 6c tan"' 6w 



30(27r)'^mfjm^ 



5m[[5mY — m^) 



2\5/2 



2 g^Gj, tan' 9 



w 



6m{{6mf - m^Y^\2{6mf - 23m^), 



3 5'^(27r)^m^^m^ 



(7) 
(8) 
(9) 



where the notations and the relevant Lagrangian, the exact decay rates are given in Appendix 
A. Notice that when 6m > mr the two-body decay is open and the lifetime of the stau, Tf^, is 
< 10~^^ sec. However, for 5m < m^ this decay is closed and the three or four-body decays 
are suppressed at least by an additional ((5m)^G|.(/^/mT-)^l/(30(27r)^) ~ 10~^^ with 6m ~ 2 
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GeV. Therefore the stau becomes long-hved and the phenomenology of the MSSM changes 
dramatically. [63] 



III. LFV AND LONG-LIVED SLEPTON 

In the previous section we have seen that 5m < rrir is indeed possible in the framework 
of a CMSSM without LFV couplings. The lifetime of the NLSP stau in this scenario is 
increased by many orders of magnitude. However, a complete absence of flavour mixing in 
the sfermion mass matrices is not expected in realistic models. In general, we would expect 
a certain degree of intergenerational mixing to be present in the slepton sector. Once a 
new source of LFVs is introduced, the NLSP two-body decay channels into electron and/or 
muon can open again. In this case, the lifetime is inversely proportional to the square of 
the mixing of selectron and smuon with stau and therefore the measurement of the lifetime 
shows a strong sensitivity to LFV parameters. In this section, we show the dependence of 
the lifetime of the lightest slepton on the right-handed and the left-handed slepton mixings. 
All the numerical results of lifetimes and figures presented below are calculated using the 
exact formula in Appendix A. 

To understand the dependence of lifetimes on LFV parameters, it is convenient to intro- 
duce the so-called Mass Insertions (MI), 11)0/3, defined 

where a, P denote the lepton flavours. M^y^^^ and M^^j^^ are diagonal elements of the 
slepton mass matrix, and AM'^^jj^^^^ are off-diagonal elements we introduced. In terms of 
these mass insertions, the two-body decay rate is approximately given by 

r2-body = ^{5rnn\gt,? + \9U\ (H) 

where a = e, fi. can be approximated in the mass insertion as shown in Fig. 3. In the 
case of right-handed slepton mixing, we have, 

/If e T\4e 

^^L^O, ^^L^tan^H^— (12) 

while in the left-handed slepton mixing case, these couplings are given 

r ^l,n m.(Ao-/itan/j) M^.M^, „ ^ 

To analyze the effects of the presence of a non-vanishing leptonic mass insertion on the 
NLSP lifetime we choose four points with different mass differences as shown in Table II. 
In Fig. 4 (a), we show the lifetime of the lightest slepton, tj^, as a function of (5|jii;)er, 
encoding the right-handed selectron-stau mixing, that we vary from 10~^° to 10~^. We can 
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FIG. 3: MI Feynman diagrams of two-body slepton decay: The diagram on the left side depicts 
two-body decay in the presence of and the diagram on the right side in the presence of S'lj^. 
The circle-crosses represent mass insertions for flavor and left-right chirality changes. Propagators 
for intermediate states are shown below corresponding scalar lines. 



No. 


5m (GeV) 


m^o (GeV) 


mj^ (GeV) 




(X 10-10) 


A 


2.227 


323.1549 


325.3817 


0.110 


10.32 


B 


1.650 


325.5601 


326.2147 


0.102 


10.25 


C 


0.407 


327.6294 


328.0365 


0.085 


10.09 


D 


0.092 


328.4060 


328.4981 


0.081 


10.06 



TABLE II: Table of the mass difference and the lightest slepton, neutralino masses, mo, and 
tan/3 are fixed to 260 GeV, 600 GeV and 30, respectively. The values of neutralino abundance and 
are shown for the reference. 



see that the hfetime for Sm > (case A in Table. II) does not change, because the decay 
of slepton to tau and neutralino is always the dominant decay mode and the lifetime is 
insensitive to (< 10"^) . On the other hand, for 6m < m,-, the lifetime grows more 
than 13 orders of magnitude in the limit (5|j^)er 0, where the three- or four-body decay 
processes are dominant. Then, the two-body decay into r and Xi is forbidden but those into 
e (or /i for {Sf^p()fj^r 7^ 0) and Xi ^-^e allowed through LFV couplings. The lifetime decreases 
proportionally to |(5|j.^)eT|~^ when the two-body decay dominates total decay width. For 
values of 5^ > 10"^ the mass of the lightest slepton would be substantially changed by 
this large off-diagonal entry and this would reduce the mass difference changing the simple 
|((5|j.^)er|~^ proportionality. This can be seen as an increase of the lifetime in the case D at 
(5|j.^)^T- ^ 10~^ in Fig. 4 (b). Although S'^j^ > 10"^ is still allowed by experiments as shown 
in Table III, we concentrate our discussion on Sf/^j^ < 10"^ to show the sensitivity of this 
process to small S'^pfS. In fact, from these figures we can see that the lifetime indeed has 
very good sensitivity to small 5^s in this scenario. The present bounds on the different mass 
insertions for the spectrum considered here are shown in Table III. Comparing this table 
with the figures we can see that the sensitivity obtained in the measurement of the NLSP 
lifetime can not be reached by indirect experiments as r — >^ /i 7 and r — > e 7, etc. Plot (b) in 
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FIG. 4: The lifetime of the hghtest slepton as a function of S'^j^. The left panel, (a), is the lifetime of 
the lightest slepton with the right-handed selection and stau mixing, and the right panel, (b) , is the 
one with the right-handed smuon and stau mixing. In both panel, the solid (red), dashed (green), 
dotted (blue), dotted-dashed (pink) line correspond to 6m = 2.23 (A), 1.65 (B), 0.41(C), 0.09 {D) 
GeV, respectively, mo = 260 GeV, = 600 GeV and M1/2 is varied. 0.08 < n^oh^ < 0.14 is 
required. These lifetimes are calculated with the exact formulas in the Appendix A 
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0.33 


0.21 
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0.00030 
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0.043 


RR 


0.0060 
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0.0019 


0.44 
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0.0012 
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TABLE III: Mass insertion bounds for a typical point in the long-lived stau coannihilation region 
with mo = 267 GeV, M1/2 = 750 GeV, ^0 = 600 GeV and different values of tan/3 [48, 50]. These 
bounds are obtained from the experimental limits {BR{fi — > 67) < 1.2 x 10~^^, BR(t 67) < 
1.1 X IQ-\BR{t ^7) < 0.45 X 10"^). 



Fig. 4 corresponds to the right-handed smuon-stau mixing, (5|j/j)^t- The lifetime is constant 
for 5m = 0.09 GeV in small {S'^fi)ij.T as the two-body decay in muon is still forbidden. For 
larger dm, we observe again the same dependence on (5^i^)/iT- 

Figs. 5 (a) and (b) show the lifetime of the lightest slepton as a function of (5£i)eT and 
(5£^)^^, respectively. The different curves correspond to the same mass differences used in 
Fig. 4. The dependence on these mass insertions in these figures is completely analogous 
to the behavior observed in Fig. 4. In fact, both figures would be identical if we replace 
(^ll) by 10 X (^R/j). This is due to the fact that in this case, the lightest slepton decays 
into the electron or the muon through the left-handed stau component. The mixing of right 
and left-handed staus, Eq. (13), is proportional to mr{Ao — /itan/9)/(M|.^^ — Ml^). In the 
region of parameter space we are considering this factor is approximately 0.1. Thus, we 
can see that, also in the case of (SIl), the lifetime is sensitive to the presence of very small 
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FIG. 5: The lifetime of the hghtest slepton as a function of The Unes and the parameters are 
the same as Fig. 4. 



lepton flavour violating couplings and therefore to the presence of mass insertions orders of 
magnitude smaller than the present bounds. 

At this point, we have seen that even very small values of the mass insertions are 
observable if nature happens to choose the coannihilation region in the MSSM for rela- 
tively large values of Mi/2- Now, we can ask what is the expected size for these 5^ in 
different models. The first example, always present if we introduce right-handed neutri- 
nos in the theory, is the supersymmetric seesaw mechanism. In this model we simply in- 
troduce right-handed neutrinos with Yukawa couplings to the Higgs and lepton doublets 
[31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45]. These Yukawa couplings induce 
some non-vanishing off-diagonal entries in the left-handed slepton mass matrices through the 
renormalization group evolution (RGE) from Mqut to My/- The value of the induced S'^s 
depends on the considered neutrino Yukawa couplings and this is model dependent. Never- 
theless, it is usually assumed that the neutrino Yukawa couplings are related to the up-quark 
Yukawa couplings, as happens in 5*0 (10) GUT models. Then, we have two main possibilities 
for this Yukawa matrix, it can have large (MNS-like) or small (CKM-like) mixings [40]. In 
the case of large mixings with 1/3 ~ i/^ ~ 1, where and yt are the largest eigenvalues in 
the neutrino and up-squark Yukawa matrices, (5£i)eT — 0.1 t/es, with U^z < 0.18 still not 
measured in the MNS matrix, and (5£i)/xr — 0.05. This mixing is very large and would 
predict a lifetime < 10~^^ sec. Even in the case of small mixings (CKM-like), the induced 
deltas are still sizable, (5£i)er — 0.0008 and (5£i)^T — 0.004, corresponding to lifetimes, 

~ 10~^° to 10~^^ sec and r^^ ~ 10~^^ to 10~^^ sec respectively. These RGE effects from 
the seesaw mechanism are always present in the model irrespective of the initial structure 
of the slepton mass matrices at the GUT scale. It is also possible that the slepton mass 
matrices have a non-trivial structure at Mqut if the same mechanism that generates the 
flavour structure in the Yukawa couplings generates simultaneously a flavour structure in 
the soft-terms. This is the case of flavour symmetries in SUSY [11, 46, 47, 48, 49]. As an ex- 
ample we take the 5'f/(3) flavour symmetry [46, 48]. Typical values of the deltas in this case 
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are ((5|j^)er ^ 0.001, (5|j^)^r ^ 0.02, (5ii)er ^ 0.0008 and ((5ii)^r ^ 0.004 and therefore 
expected lifetimes would be in the range 10~^° to 10"^^ sec. Notice that this model predicts 
only moderate flavour changing couplings in the slepton mass matrices and values in other 
models are typically larger or the same order as this [11]. Non- vanishing off-diagonal entries 
in the slepton mass matrices are also generated by SU{5) or 5*0(10) RGE evolution from 
Mpianck to Mqut [51, 52, 53, 54], and similar sizes of d'^s are expected in these models. 

IV. LHC PHENOMENOLOGY 

In this section, we discuss the expected phenomenology at LHC experiments, focusing 
mainly on the ATLAS detector [55, 56], of the long-lived slepton scenario. The lightest 
slepton is the NLSP, and therefore a large number of sleptons is expected to be produced 
via cascade decays of heavier SUSY particles. When the long-lived sleptons have lifetimes 
between 10~^ and 10~^^ sec, we will have a chance to observe decays of the slepton inside 
the detector. The ATLAS detector, whose dimensions are 25 m in hight and 44 m in length, 
consists of inner detector, calorimeter and muon detector [55] . The inner detector surrounds 
the collision point and is contained within a cylindrical envelop of length ±3.51 m and of 
radius 1.15 m. In the inner detector, the pixel detector, the semiconductor tracker (SCT) and 
the transition radiation tracker (TRT) are installed to measure momenta and decay vertices 
of charged-particles. The pixel detector is placed around the collision point and consists of 
three barrels at average radii of 5 cm, 9 cm and 12 cm, and three disks on each side, between 
radii of 9 and 15 cm. The SCT surrounds the pixel detector and consists of four barrels at 
radii of 30, 37, 44 and 52 cm. The TRT is installed on outside of the SCT. It covers radial 
range from 56 to 107 cm. Surrounding the inner detector, hadronic and electromagnetic 
calorimeters are placed over 4 m radius and 8.4 m length. The calorimeters are surrounded 
by the muon spectrometer which covers the ATLAS detector. In the following, we show 
that, in the relevant region of SUSY parameter space, we can see decays of sleptons with 
the lifetime up to 10~^ seconds within the ATLAS detector. Although we discuss only the 
ATLAS detector, analysis is similar to the CMS detector [57]. 

First, we estimate the number of lightest sleptons produced at LHC. In our scenario, 
the lightest slepton is mainly a right-handed stau and the lightest neutralino is mainly 
bino. Then, the right-handed staus are mainly produced via decays of first two generation 
left-handed squarks, (number in parenthesis corresponds to branching ratios in each decay 
[58, 59]) 

qL^xt + q (0.63), (14) 
Xt^h + ^r (0.64), p, + T (0.27), (15) 
v^^l-^ + W (0.82), (16) 

qL^xl + q (0.36), (17) 
Xl^h + r (0.66), Vr + Vr (0.25), (18) 

where g, xf and x% denote the SM quarks, the lighter charginos and the second lightest neu- 
tralinos and z/^, W and z/^ denote tau neutrino, weak boson and tau sneutrino, respectively. 
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The 3rd generation squarks have many different decay chains. For example 





-^Xf + b, (0.20, 0.25) 


(19) 




Xt h + Z^r, 


(20) 


tl,2 


^Xl + t, (0.10, 0.10) 


(21) 




xl h + 


(22) 


&1,2 


^Xt + t, (0.36, 0.12) 


(23) 






(24) 




^X2 + ^ (0.20, 0.10) 


(25) 




Xl h + l^r- 


(26) 



Total branching ratios for li production are 0.86 in decay, 0.72/0.90 for ti/t2 and 0.87/0.67 
for 61/62 at tan /3 = 35 [58, 59, 60]. On the other hand, the right-handed squarks, couple only 
to bino and higssino, having small Yukawa couplings, decay almost completely to the lightest 
neutralino and quarks. Therefore branching ratios of them into li and quarks/leptons are 
negligible. The total cross section of SUSY pair production in our scenario is given in [61], 

(TsusY = 130 fb. (27) 

When we assume the integrated luminosity, C-mt = 30 fb~^, which corresponds to three year 
of data taking with the luminosity, C = 10^^ cm~^ s~^, the total number of SUSY pairs 
is 3900. Then, squarks are pair-produced with equal probability for left and right-handed 
squarks, thus, the number of lightest sleptons produced, Nj_^ , is estimated as 

Nj^ ^ 4290. (28) 

Therefore we would have 4290 long-lived sleptons produced in ATLAS that could decay 
inside the detector depending on the lifetime. The decay probability is given 

Pdec(0 = l-expf-— (29) 

where j3 = v/c and 7 = 1/ ^Jl — [3"^ with c the speed of the light, / is the distance between 
production and decay point. To obtain the expected number of decays, we need to know the 
distribution of /57. This would require a full analysis with Monte-Carlo simulation and it is 
beyond the scope of this paper. For the spectrum shown in Table I, the typical momentum 
of the lightest slepton is expected to be between 500 and 900 GeV. /37 corresponding to this 
range of momenta is 

1-53 < /37 < 2.75, (30) 

hence we can assume /37 to be 2. Then, using Eq.(29), the expected number of the decays, 
A'^dec, within a distance / is given 

N^^^il) = A^ Pdec(/) = A^^, (1 - exp (-;^) ) • (31) 
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175.1 


10"^ sec. 


3.6 
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1731.0 


2658.3 


4223.5 


10-9 sec. 


343.0 


2425.6 


4265.5 


4289.7 


4290.0 


10-1° sec. 


2425.6 


4289.0 


4290.0 


4290.0 


4290.0 


10-11 sec. 


4289.0 


4290.0 


4290.0 


4290.0 


4290.0 


10-12 sec. 


4290.0 


4290.0 


4290.0 


4290.0 


4290.0 



TABLE IV: The expected number of slepton decay in the ATLAS detector. The length is the 
minimum (5 cm) and medium (50 cm) distance to the pixel detector and the maximum distance to 
the outer boundary of the detectors from the interaction point and corresponds to pixel detector 
(3.1 m), calorimeter (5.8 m) and muon spectrometer (25.0 ni). The lifetime is varied from lO^i^ 
to 10"^ seconds. The number of sleptons produced at LHC assuming the integrated luminosity, 
£int = 30 fb-i is 4290 and is fixed to 2. 

In Table IV, we show the expected number of slepton decays in each detector, assuming 
the integrated luminosity, Ci^t = 30 fb-i and P'j = 2. As can be seen, when the lifetime, r^^, 
is below IQ-^ sec, most of the sleptons decay inside the pixel detector. When r^-^ ~ IQ-^ 
sec, almost half of them decay inside the inner detector and nearly all of them decay within 
the detector. If r^^ is between IQ-*^ and 10"^ sec, several hundreds of slepton decays occur 
inside the ATLAS detector. Almost all of them escape from the detector when r^^ > IQ-^ 
sec, although we expect of the order of 10 decays inside the detector for tj^ ~ IQ-^ sec. 

Sleptons with different lifetimes would give different signatures in the ATLAS detector. 
Depending on the lifetimes, sleptons would decay inside the detector or escape the detector 
as a stable heavy charged-particle. For r^^ < lO^n sec, since almost all of the sleptons would 
decay before they reach the first layer of the pixel detector, we would not observe any heavy 
charged-particle track in the detector. In this case, it would be more difficult to identify the 
presence of long-lived sleptons at the ATLAS detector. This problem will be addressed in 
a future work [62]. For tj_^ ~ IQ-i^ to IQ-^ sec, almost all the sleptons would decay inside 
the pixel detector and leave a charged track with a kink. Then, a different charged-particle 
would cross the outer detectors and a corresponding track and/or hit would be seen in each 
detector. Thus, by combining with signals in the outer detectors, we could identify whether 
the outgoing charged-particle is an electron or muon. Then, if we can fix the mass difference 
between the lightest slepton and neutralino, we can determine the value of the mass insertion 
parameters from Figs. 4 and 5. In the case of right-handed slepton mixing case, lifetimes 
between IQ-i^ and IQ-^ sec. would correspond to (^|j/{)eT between IQ-'' and IQ-^ with the 
mass difference, rrie < 5m < rrir, and {S'^fi)fiT between IQ-^ and IQ-^ with < 6m < mr- 
Similarly, in the case of left-handed slepton mixing, the same lifetimes would correspond to 
i^LL)er bctwceu 4 X IQ-^ and IQ-'^ with m^ < Sm < m^-, and (5|,i)/iT between 4 x IQ-^ and 
IQ-^ with m^ < 6m < m^. For r^^ ~ IQ-^ sec, about half of sleptons would decay inside the 
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inner detector and the rest would decay inside the calorimeters and/or muon spectrometer. 
In this case, it would be important to determine whether the decay is LFV two-body decay 
or lepton flavour conserving three-body decay [62]. For lifetimes between 10"'' and 10~^ sec, 
very few sleptons would decay inside the pixel detector and most of them would escape the 
detector leaving a charged track with a corresponding hit in muon spectrometer. Even for 
particles escaping the detector, we could use muon spectrometer to determine slepton mass 
and momentum as studied in [21, 22]. For r^^ > 10"^ sec, very few sleptons would decay 
inside the ATLAS detector. In this case, we can put lower bound on the lifetime. Then, in 
a gravity-mediated scenario, this would correspond to a stringent upper bound on (5^s. 

This scenario of long-lived slepton is in principle similar to the situation in gauge- mediated 
SUSY breaking scenario [23]. For similar lifetimes, more detailed analysis is needed to 
distinguish scenarios. However most of sleptons in gauge-mediated SUSY breaking scenario 
would escape from the detector and we could see only a few events. 

V. SUMMARY 

In this work, we have studied lepton flavour violation in a long-lived slepton scenario 
where the mass difference between the NLSP, the lightest slepton, and the LSP, the lightest 
neutralino, is smaller than the tau mass. With this small mass difference, the lifetime of 
the lightest slepton has a very good sensitivity to small lepton flavour violation parameters 
due to the more than 13 orders of magnitude suppression of the three-body decay rate with 
respect to the two-body decay rate. 

We have shown that this small 5m is possible even in the framework of the CMSSM. 
We selected the region in the CMSSM parameter space with 5m < mr, taking into account 
correct dark matter abundance and direct bounds on SUSY and Higgs masses and b s 'j 
constraint. We have used the experimental results of the anomalous magnetic moment of the 
muon to indicate the favored regions at given confldence level. We have found that a small 
5m, region consistent with all experimental constraints lies on sizable part of the parameter 
space at M1/2 ^ 700 GeV, m^ = 100 to 600 GeV for different values of tan/3. We have 
also shown that the most favored region, which is consistent with at 2 cr, corresponds to 
240 < mo < 260 GeV and 700 < M1/2 < 800 GeV at tan/3 = 35. 

Then, we have analyzed the dependence of the lightest slepton lifetimes on different mass 
insertions, ('5^^/x,x,)eT, ^1T for values of 5m from 2.23 to 0.09 GeV. We found that the lifetimes 
are proportional to l^*^]"^ until three- or four-body decays become comparable. There is a 
difference of approximately a factor of 10 in the sensitivity of the lifetime on and 511^. 
This is due to the different proportion of left-handed and right-handed staus in the lightest 
slepton. By comparing the values of the bounds on ^'^s shown in Table III, we can see that, 
in this scenario, the lifetimes are sensitive to much smaller values of these (5^s, even to future 
sensitivities of proposed experiments. 

Finally, we have discussed the expected phenomenology at LHC experiments, mainly 
concentrating on the ATLAS detector. We have estimated the number of slepton decays 
in the different detectors, assuming an integrated luminosity Ci^t = 30 fb~^ and iS'j = 2 
(Table. IV). We have seen that the ATLAS detector can observe lifetimes in the range of 
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10~^^ to 10~^ sec, and these lifetime would correspond to ((5^R)eT,Atr between 10~^ and 10~^ 
and {51i)eT,iiT between 4 x 10~^ and 10~^. Therefore we have shown that in the long-lived 
slepton scenario, the LHC offers a very good opportunity to study lepton flavour violation. 
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APPENDIX A: INTERACTION LAGRANGIAN 

In this appendix, we give the relevant Lagrangian and the decay widths in the presence 
of LEV. 

The slepton-neutralino-lepton interaction Lagrangian in neutralino mass eigenstate basis 
and sfermion flavour basis is 

C^_^o_i = V2g, ^ [{gt'^J: + gt'^t^MPLla + (^r/:+3 + 9^^~foMPRla] + h.c. (Al) 

where denote mass eigenstate neutralinos and i runs over 1 to 4, and fa = (e^, /il, tl) and 
/a+3 = (e_R, yUR, tr) denote flavour sleptons and a runs over 1 to 3. Pl,r are the chirality 
projection operators and g2 is the SU(2) coupling constant. The constants, gf^, couple 
sleptons of chirality "a" to leptons of chirality "6". 

gf^ = -^*[T,l{N),, + (g - T,l) i^riOwiNU. 9^' = V* J Nh, (A2) 

niw cos p 

^f/ = r/,gtan^^(iV),i, ^fj^ = (N^, (A3) 

mw cos p 

where Q and T^^ is the electric charge and SU{2) charge of leptons, and 9w is the Weinberg 
angle. Nij is the diagonalization matrix of the neutralino mass matrix. 
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The slepton mass matrix is diagonalized by a unitary matrix and the mass eigenstate 
sleptons, Ik, with masses are nij^, are expressed by 

h = Rkafa- (A4) 

We define R according to mj^ > mj^^ for k > k'. Then, the relevant interaction Lagrangian 
in the mass eigenstate basis is given by 

4 6 
i=l k=l a=e,fi,T 

G AG 
+ -^Y^ i^ali.PLQJ'' + -1=- Yl {^o.li^Qi.^Pli^h) + h.c. (A5) 

* a=e,fi,T ^ a,f5=e,fi,T 

Here Gi? is the Fermi couphng constant, is the pion current, and g^^^, are 

9iak = 9iaRka + g^^Rka+S, (A6) 
9iLk = 9iL^Rka+3 + gi^Rka- (A7) 

The expression of these couphng constants in MI approximation is given in Eqs. (12) and 
(13). 

The decay widths can be obtained by replacing g^ and gn in Ref. [2] with gi^i- The 
two-body decay width, /i ^ Xi + where la = (e, /i), is given by 

r2-body = gj^3 lsm{6m + 2m^o){5m{6m + 2m^o) - 2mfJ + mf] 

h 

><{i\9tif + \9ti\')iSmi5m + 2m^o)-ml)-ARe[gti9L*]'^^^ 

(A8) 

The three-body decay width, h ^ Xi + + n where la = (e, fi, r), is given by 



_ glGlP^ cos' 9, {{6m)^ -ml) 
8(27r)3m? 



(qz — ml 



2\2 



X 



^ {\9iaihl + \9iai\^rnl){{6m)' + 2m^o6m - ql) - Re[g^^^gti]m^omi^ql 

(A9) 

where = {5mY — ((5m)^ — m^)x, and F; is the decay width of the charged lepton, / 
(Fe = 0). 6c is Cabbibo angle. 

The 4-body decay width, h Xi + ^la + ^ifs + where la = (e, /i, r) and Ip = (e, /i), 
is given by 



4-body 



giGl [{6my-ml) 
' 12(27r)5m? 
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X 



^(l^iail Qif, + l^iail ^ijii^^) + 2m^o5m - g^^) - Re[^i^i^i^i ]m^om,^g;^ 



X < 12mf oMog 



(AlO) 



where = (^m)^ — ((^m)^ — m^^ )x. 
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